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Direct imaging of exoplanetary systems is a powerful technique that can re- 
veal Jupiter-like planets in wide orbits, can enable detailed characterization of 
planetary atmospheres, and is a key step towards imaging Earth-like planets. 
Imaging detections are challenging due to the combined effect of small angu- 
lar separation and large luminosity contrast between a planet and its host star. 
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High-contrast observations with the Keck and Gemini telescopes have revealed 
three planets orbiting the star HR 8799, with projected separations of 24, 38, 
and 68 astronomical units. Multi-epoch data show counter-clockwise orbital 
motion for all three imaged planets. The low luminosity of the companions 
and the estimated age of the system imply planetary masses between 5 and 13 
times that of Jupiter. This system resembles a scaled-up version of the outer 
portion of our Solar System. 

During the past decade various planet detection techniques - precision radial velocities, 
transits, and microlensing - have been used to detect a diverse population of exoplanets. How- 
ever, these methods have two limitations. First, the existence of a planet is inferred through 
its influence on the star about which it orbits; the planet is not directly discerned [photomet- 
ric signals from some of the closest-in giant planets have been detected by careful analysis of 
the variations in the integrated brightness of the system as the planet orbits its star [1]]. Sec- 
ond, these techniques are limited to small (transits) to moderate (precision radial velocity and 
microlensing) planet-star separation. The effective sensitivities of the latter two techniques di- 
minish rapidly at semi-major axes beyond about 5 AU. Direct observations allow discovery of 
planets in wider orbits and allow the spectroscopic and photometric characterization of their 
complex atmospheres to derive their physical characteristics. 

There is indirect evidence for planets in orbits beyond 5 AU from their star. Some images of 
dusty debris disks orbiting main- sequence stars (the Vega phenomenon) show spatial structure 
on a scale of tens to hundreds of astronomical units [2]. The most likely explanation of such 
structure is gravitational perturbations by planets with semi-major axes comparable to the radius 
of the dusty disks and rings [see references in jSl]. 

The only technique currently available to detect planets with semi-major axes greater than 
about 5 AU in a reasonable amount of time is infrared imaging of young, nearby stars. The de- 
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tected near- infrared radiation is escaped internal heat energy from the recently formed planets. 
During the past decade, hundreds of young stars with ages < 100 Ma have been identified within 
~100 pc of Earth [|4l|5l, and many of these have been imaged in the near-IR with ground-based 
adaptive optics (AO) systems and with the Hubble Space Telescope. Direct imaging searches 
for companions of these stars have detected some objects that are generally considered to be 
near or above the mass threshold 13.6 Mjup dividing planets from brown dwarfs [see [16| for an 
example and Q for a list of known substellar objects orbiting stars], and one planetary mass 
companion that is orbiting a brown dwarf, not a star [[8]|. Recently, Lafreniere et al. flU have 
detected a candidate planet near a young (5 Ma) star of the Upper Scorpius association, but a 
proper motion analysis is required to confirm that it is bound to the host star and not an unre- 
lated low-mass member of the young association. In this issue, Kalas et al. report the detection, 
in visible light, of a candidate planetary mass companion near the inner edge of the Fomalhaut 
debris disk ifTOl . Non-detections of the candidate companion at near-IR wavelengths suggest 
that the detected visible flux may be primarily host-star light scattering off circumplanetary dust 
rather than photons from the underlying object. A statistical Bayesian analysis of a dedicated 
AO survey of nearby young F-, G- and K-type stars shows that exoplanets are relatively rare at 
separations > 20 AU around stars with masses similar to the Sun [fTTl . 

Bright A-type stars have been mostly neglected in imaging surveys since the higher stel- 
lar luminosity offers a less favorable planet-to-star contrast. However, main sequence A-type 
stars do have some advantages. The higher-mass A stars can retain heavier and more extended 
disks and thus might form massive planets at wide separations, making their planets easier to 
detect. Millimeter interferometric continuum observations of the nearest Herbig Ae stars, the 
precursors to A-type stars, indicate that these are encircled by disks with masses up to several 
times the Minimum Mass Solar Nebula [12J, the minimum amount of solar abundance material 
(O.OlMsun) required to form all planets in the Solar System |[T3l . Associated millimeter line 
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observations resolve these gas disks and indicate that their outer radii are 85-450 AU [[T2l . The 
most exceptional example of a young A-star disk is the one orbiting IRAS 18059-321 1, which 
is estimated to have a mass of 90 times the Minimum Mass Solar Nebula and an outer radius 
extending to ~3,000 AU [fT4ll . Radial velocity surveys of evolved A stars do seem to confirm 
these hypotheses by showing a trend of a higher frequency of planets at wider separations ifTSl . 
In this article, we describe the detection of three faint objects at 0.63", 0.95" and 1.73" (24, 38 
and 68 AU projected separation, see Fig. 1) from the dusty and young A-type main sequence 
star HR 8799, show that all objects are co-moving with HR 8799, and describe their orbital 
motion and physical characteristics. 

HR 8799 Stellar Properties 

HR 8799 (also V342 Peg & HIPl 14189, located 39.4 pc US from Earth) is the only star known 
that has simultaneously been classified as 7 Doradus (variable), A Bootis (metal-poor Population 
I A-type star) and Vega-like (far-IR excess emission from circumstellar dust) ifTTlfTSl . A fit to 
the Infrared Astronomical Satellite (IRAS) and Infrared Space Observatory (ISO) photometry 
indicates that it has a dominant dust disk with temperature of 50 K |l3l [TH . Such blackbody 
grains, in an optically thin disk, would reside ~75 AU (~ 2") from HR 8799. This would place 
the dust just outside the orbit of the most distant companion seen in our images (see Fig. 1), 
similar to the way the Kuiper Belt is confined by Neptune in our Solar system. 

The fractional IR luminosity (LiR/Lboi = 2.3X10""^) EOl |T9ll is too bright to come from a 
geometrically thin, flat disk orbiting at such large distances from HR 8799. Such an optically 
thin disk would need to be warped or puffed up in the vertical direction, plausibly by the grav- 
itational influence of nearby planets. Submillimeter photometry indicates a dust mass of 0.1 
Earth masses (EH, making it one of the most massive debris disks detected by IRAS llT9l . 

When planets form, gravitational potential energy is released and turned into heat in their 
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interior. As planets do not possess any internal nuclear energy source to maintain their tem- 
perature, they cool down and become less luminous with time. For massive planets, this self- 
luminosity can dominate over their stellar insolation for hundreds of millions or billions of 
years. With some assumptions on the initial conditions at the time of formation, a planet's mass 
can be derived simply by estimating the planet's luminosity and the system age. Our age esti- 
mate for HR8799 is based on four lines of evidence: the star's galactic space motion, the star's 
position in a color-magnitude diagram, the typical age of A Boo and 7 Dor class stars and the 
large mass of the HR 8799 debris disk. 

Most young stars in the solar neighborhood have Galactic space motions (UVW) that fall in 
limited ranges. HR 8799's space motion with respect to the Sun, as calculated from published 
distance, radial velocity and proper motion, is UVW = (-1 1.9, -21.0, -6.8 km s~^) [[T6ll22l . This 
UVW is similar to that of other stars with an age between that of the TW Hydra association 
[8 Ma a [51] and the Pleiades [125 ± 8 Ma [lH]. The UVW of HR 8799 is similar to that 
of members of the 30 Ma old, southern hemisphere, Columba and Carina Associations |l5l. 
Calculations of the UVW of the young stars HD 984 and HD 221318, which lie near HR8799, 
show that their space motions are similar to that of HR8799. We estimate the ages of HD 984 
and HD 221318 to be 30 and 100 Ma, respectively, while the FEPS team estimates the age of 
HD 984 to be 40 Ma EH- Overall, the UVW of HR 8799 is clearly consistent with those of 
young clusters and associations in the solar neighborhood. Of course, in this UVW range of 
young stars, there are also older stars with random motions; so other, independent, methods 
must also be employed to place limits/constraints on the age of HR 8799. 

HR 8799 is also found below the main sequence of the Pleiades, a Per (70 Ma) and IC2391 
(50 Ma) on a Hertzsprung-Russell diagram. This is consistent with a younger age compared 
to that of the Pleiades [|25l . Even with the more recent Tycho measurement and correcting for 
the star's low metallicity, so that B-V is increased and lies between 0.26 and 0.3, HR 8799 
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still lies low on the Hertzsprung-Russell diagram when plotted against known young stars ll25l . 
consistent with our young age estimate. 

The A Boo stars are generally thought to be young, up to a few 100 Ma ll26l . The 7 Dor class 
stars are probably also young; they are seen in the Pleiades and in NGC 2516 (age ~100 Ma), 
but not in the Hyades (age ~650 Ma) [l27ll . 

Finally, the probability that a star has a massive debris disk like HR 8799 declines with age 
[fT9ll . Considering all of the above, we arrive at an estimate of 60 Ma and a range between 
[30,160] Ma, consistent with an earlier independent estimate of 20 — 150 Ma [20J. The con- 
servative age upper limit for HR 8799 is chosen to be the ~ 5a upper limit to the Pleiades 
age. 

Observations 

The sensitivity of high-contrast ground-based AO imaging is limited primarily by quasi-static 
speckle artifacts; at large separations (> 0.5"), the main source of speckles is surface errors 
on the telescope primary mirror and internal optics. To remove this noise, we used angular 
differential imaging (ADI) in our observations [|28ll29l . This technique uses the intrinsic field- 
of-view rotation of altitude/azimuth telescopes to decouple exoplanets from optical artifacts. 
An ADI sequence is obtained by keeping the telescope pupil fixed on the science camera and 
allowing the field-of-view to slowly rotate with time around the star. Our observations were 
obtained in the near-infrared (1.1 to 4.2 microns), a regime where the planets are expected to be 
bright and where the AO system provides excellent image correction. ADI sequences at various 
wavelengths were acquired using the adaptive optics system at the Keck and Gemini telescopes 
and the corresponding facility near-infrared cameras, NIRC2 and NIRI, between 2007 and 2008 
(see Fig. 1). During each observing sequence, we typically obtained a mix of unsaturated short- 
exposure images of the star, to determine its precise location and brightness, together with a 
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set of 30-s exposures that overexposed the central star but had maximum sensitivity to faint 
companions. Some coronagraphic images were also acquired with NIRC2 to benefit from the 
simultaneous photometric calibration achievable with a partially transmissive focal plane mask. 
The b and c companions were first seen in October 2007 Gemini data; the d component was 
first detected in Keck data in 2008. The b and c components were also visible in a re-analysis 
of non-ADI Keck data obtained for a related program in 2004 [the data sets and the reduction 
technique are described further in the supplemental online material (SOM)] 

Astrometric Analysis 

Following the initial detection of the companions, we evaluated their positions relative to the star 
to confirm that they are co-moving with it (possibly including orbital motion) and not unrelated 
back/foreground objects (see SOM; Table 1, Table S2, and Fig. 2). Because HR 8799b was 
visible in the 2004 Keck images, we have more than 4 years of time baseline for proper motion 
measurements. With the large proper motion of HR 8799 (0.13'7year), the HR 8799b object is 
shown to be bound at a significance of 98 times the estimated 1-sigma uncertainty. Additionally, 
the data show that it is orbiting counter-clockwise. It moved 25 ± 2 mas/year (0.98 AU/year) 
southeast during the 4 year period. Its detected orbital motion is near perpendicular to the line 
connecting the planet and primary, suggesting that the system is viewed nearly pole-on and 
that the orbit is not very eccentric. The near face-on perspective is further supported by the 
slow projected rotational velocity of HR 8799 [~40 km sec~^ [fTTl l; this is well below average 
for late-A and early-F type stars lf30l . If we assume that it has a semi-major axis of 68 AU, a 
circular orbit, a pole-on view, and a host stellar mass equal to 1.5 Solar masses, then the orbital 
period and motion of HR 8799b are ~450 years and 0.93 AU/year (24 mas/year) respectively, 
consistent with our measurements. 

HR 8799c is also detected, at lower significance, in the 2004 data set. The measurement of 

7 



its 4 year proper motion confirms that it is bound to the star at the 90 sigma level. Its orbit is 
also counter-clockwise at 30 ± 2 mas/year (1.18 AU/year). For its semi-major axis of 38 AU, 
the orbital period is ~ 190 years and the expected orbital motion is 1.25 AU/year (32 mas/year). 
Again, the orbital motion is close to being perpendicular to the line connecting the planet to the 
primary. 

HR 8799d was first detected in the July 2008 data set. The two months of available proper 
motion measurements are sufficient to confirm that it is bound to the star at the ~6 sigma level. 
The available data is also consistent with a counter-clockwise orbital motion of 42 ±27 mas/year 
(1.65 AU/year). For a semi-major axis of 24 AU, the orbital period is 100 years and the expected 
orbital motion is 1.57 AU/year (40 mas/year). 

HR 8799bcd Photometric Analysis 

AU three companions are intrinsically faint and have red near-IR colors that are comparable to 
those of substellar-mass objects with low effective temperatures (see Table 1). Compared to 
old field brown dwarfs (objects with masses between planets and stars), all three companions 
lie at the base of the L dwarf spectral sequence - objects known to be cool and have dusty 
clouds in their atmospheres (Fig. 3). Two candidate free-floating Pleiades brown dwarfs, with 
comparable colors and absolute K-band magnitudes to HR 8799c and d, are consistent with 
a mass of ~ 11 Mjup from evolutionary models [[3T1l . If HR 8799 is (as is likely) younger 
than the Pleiades, the c and d companions would be even less massive . HR 8799b is fainter 
than all of the known Pleiades substellar members and thus is below 11 Mjup (Fig. 3). All 
three companions stand apart from the older, more massive brown dwarfs in a color-magnitude 
diagram. The known distance to HR 8799, and photometry for each companion that covers a 
substantial fraction of the spectral energy distribution (SED), allow for a robust measurement 
of the bolometric luminosity (Lboi)- We fit a variety of synthetic SEDs (generated with the 
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PHOENIX model atmosphere code) to the observed photometry for each companion assuming 
their atmospheres were either cloud-free, very cloudy, partly cloudy (50% coverage), or radiated 
like black bodies. This fitting process is equivalent to simultaneously determining bolometric 
corrections for each band-pass for various model assumptions. Luminosities were also obtained 
using the K-band bolometric corrections for brown dwarfs [|32l . Although the different models 
produce different estimates of effective temperature, the range of Lboi for each object is small 
(see Table 1), indicating that our estimate is robust against the uncertainty in the details of the 
atmosphere and clouds (see the SOM for more details). 

The cooling of hydrogen-helium brown dwarfs and giant planets is generally well under- 
stood; however, the initial conditions associated with the formation of objects from collapsing 
molecular clouds or core accretion inside a disk are uncertain. Consequently, theoretical cool- 
ing tracks of objects at young ages may not be reliable. Recent efforts to establish initial con- 
ditions for cooling tracks based on core-accretion models have produced young Jupiter-mass 
planets substantially fainter (< 10"^ Lsun) than predicted by traditional models [[33l . However, 
these hybrid models do not yet include a realistic treatment of the complex radiative transfer 
within the accretion- shock and thus provide only lower- limits on the luminosity at young ages. 
Warmer, more luminous planets originating from core accretion cannot be ruled out. 

Although HR 8799 is young, its upper age limit (~ 160 Ma) is near the time when the 
differences among cooling tracks with various initial conditions are not so dramatic and, given 
the uncertainties associated with all planet evolution models, standard cooling tracks are as 
reliable at these ages as other hybrid models. Fig. 4 compares the measured luminosities and 
age range for HR 8799 bed to theoretical hot start cooling tracks for a variety of masses [34J. 

The region occupied by all three companions falls below the lowest mass brown dwarf, well 
inside the planet regime. The masses derived from the luminosities, cooling tracks, and best 
age for bed are respectively 7, 10, and 10 Mjup. See table 1 for values of additional important 
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properties derived from the cooling track comparison, with uncertainties based on our current 
best age range. In the very unlikely event that the star is older than our estimated upper limit, it 
would need to be > 300 Ma for all three objects to be brown dwarfs. 

The large planet masses and orbital radii in the HR 8799 system are challenging to explain 
in the context of a core accretion scenario. A number of factors such as stellar mass ll35l . 
metallicity [136 L disk surface density ll37l . and planet migration in the disk ll38l influence the 
core accretion process. The stellar mass of HR 8799 is larger than the Sun. The star's metallicity 
is low, especially in refractory elements, but for a A Boo star this is usually attributed to the 
details of the stars accretion and atmospheric physics rather than an initial low metallicity for 
the system ll26ll . 

The exceptionally dusty debris disk around HR 8799 may indicate that the proto-planetary 
disk was massive and had a high surface density, factors conducive to planet formation. Alter- 
natively, the giant planets in the HR 8799 system may have formed rapidly from a gravitational 
instability in the early disk |l39l|40l. Some models HOl of such instabilities do favor the creation 
of massive planets (> 6 Mjup). 

As suggested by the color-magnitude diagram (Fig. 3), each companion appears to be at the 
edge of (or inside) the transition region from cloudy to cloud-free atmospheres. Current planet 
atmosphere models have difficulties fitting the color and spectrum features of these objects. The 
physical mechanism responsible for the clearing of clouds in ultra-cool atmospheres is not fully 
understood, but recent cloud models with vertical stratification have had some success at sim- 
ulating/producing photometric properties in this transition region [|4T]|. A modified PHOENIX 
atmospheric model was developed that incorporates cloud stratification. These updated models 
were found to match well-known brown dwarfs located in the cloudy/cloud-free transition re- 
gion. With the cloud stratification model, PHOENIX is capable of producing spectra that are 
consistent with the observed photometry and the bulk properties (effective temperature, radius. 
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and gravity) predicted by the cooling tracks (Fig. 5). Clearly these synthetic models do not 
reproduce all of the photometric data, but given the difficulty of cloud modeling, the agreement 
is sufficient to support the effective temperatures and radii determined from the cooling tracks. 

Conclusions 

The three co-moving companions of HR 8799 are significantly different from known field ob- 
jects of similar effective temperature; the only similar object known is the planetary mass com- 
panion to the brown dwarf 2M1207. Low luminosities of these companions and the young age 
for HR 8799 indicate that they have planetary masses and are not brown dwarfs. The nature of 
the system provides an additional indirect line of evidence for planetary-mass companions (and 
hence a low age). There are no known systems where multiple brown dwarfs independently 
orbit an early-type star; the only systems we know of with multiple companions in independent 
orbits are the exoplanetary systems discovered from the precision radial velocity method. Inter- 
estingly, our observations show that the HR 8799 planets orbit in the same direction, similar to 
the planets in our own solar system and consistent with models of planet formation in a disk. In 
many ways this resembles a scaled-up version of our solar system. HR 8799 has a luminosity 
of 4.9 Lsurij so the radius corresponding to a given equilibrium temperature is 2.2 times larger 
than the corresponding radius in our solar system. Because formation processes will be affected 
by luminosity - e.g. the location of the snow line where water can condense on rocky material 
to potentially form giant planet cores - one can view the three planetary companions as having 
temperature-equivalent projected orbital separations of 1 1 , 17 and 3 1 AU, to be compared with 
9.5, 19, and 30 AU for Saturn, Uranus, and Neptune. The HR 8799 planets are also consistent 
with formation through instabilities in a massive protoplanetary disk, which may form objects 
with masses above 5 Mjup f[40ll . but the core accretion scenario cannot yet be ruled out. 

The presence of these massive planets still leaves dynamic room for other Jovian-mass plan- 
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ets or even lower mass terrestrial planets in the inner part of the system. In our survey, we only 
observed a few early-type stars before making this detection, compared to similar imaging sur- 
veys of young G-, K-and M-type stars that have covered more than a few hundred targets. This 
may indicate that Jovian-mass planetary companions to early-type stars are much more common 
at separations beyond ~20 AU, consistent with what was suggested by radial velocity surveys 
of evolved A-type stars [15J. 
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Table 1: HR 8799 planetary system data 



HR 8799 


Spectral type 




A5V 




Mass 




1.5±0.3Msun 


m 


Luminosity 




4.92 ± 0.41 Lsun 


m 


Distance 


39.4 ± 1.0 pc (128 ± 3 ly) [fT6ll 


Proper motion [E,N] 


[107.88 


±0.76,-50.08 ±0.63] mas/year Ml 


Age 




60[30 - 160] Ma 




ivietaiiicity 


log([M/H]/[M/H]sun) = - 


0.47 na 


T TT T7' „ T ' 

J, H, Ks, L 


5.383 ± 0.027, 5.280 ± 0.018, 5.240 ± 0.018, 5.220 ± 0.018 




Separation w.r.t the host star in [E, N]" 




riR /yyb 


c 


a 


ZUU4 July 14 (±U.UUo ) 


[1/171 n ss/il 
[1.4 1 i, U.oo4J 


r n 7Qn n fii ol 
[— U. 1 oy, U.DizJ 




ZUU/ UCt. ZJ (±U.UUo ) 


[1.512,0.805] 


[-0.674,0.681] 




zUUo July 11 (±U.UU4 ) 


[1.527,0.799] 


[-0.658,0.701] 


[-0.208, -0.582] 


or\AO A n 10 /'_i_n nno"\ 
zUUo Aug. Iz (±U.UUz ) 


[1.527,0.801] 


[-0.657,0.706] 


[-0.216,-0.582] 


zUUo bept. lo (±U.UUo ) 


[1.528,0.798] 


[-0.657,0.706] 


[-0.216,-0.582] 


rrojectea oep. (^au j 


68 


38 


24 


uruiiai ivioiion /yearj 


0.025 ±0.002 


0.030 ±0.002 


0.042 ±0.027 


Period for pole view 








cir. orbits (years) 


~ 460 


- 190 


~ 100 


ivij (^i.z4o/im^ 


16.30 ±0.16 


14.65 ±0.17 


15.26 ±0.43 


Mh (1.633/im) 


14.87 ±0.17 


13.93 ±0.17 


13.86 ±0.22 


McH4S (1.592/im) 


15.18 ±0.17 


14.25 ±0.19 


14.03 ±0.30 


McH4L (1.681/im) 


14.89 ±0.18 


13.90 ±0.19 


14.57 ±0.23 


Mks (2.146/im) 


14.05 ± 0.08 


13.13 ±0.08 


13.11 ±0.12 


Ml' (3.776/im) 


12.66 ±0.11 


11.74 ±0.09 


11.56 ±0.16 


Luminosity (Lsun) 


-5.1 ±0.1 


-4.7 ±0.1 


-4.7 ±0.1 


Teff (K) 


870[800 - 900] 


1090[1000- 1100] 


1090[1000 - 1100] 


Radius (Rj^p) 


1.2[1.1 - 1.3] 


1.2[1.2 - 1.3] 


1.2[1.2 - 1.3] 


Mass (Mjup) 


7[5-ll] 


10[7- 13] 


10[7- 13] 
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Figure 1: HR 8799bcd discovery images after the light from the bright host star has been re- 
moved by ADI processing. (Upper left) A Keck image acquired in July 2004. (Upper right) 
Gemini discovery ADI image acquired in October 2007. Both b and c are detected at the 2 
epochs. (Bottom) A color image of the planetary system produced by combining the J-, H-, 
and Ks-band images obtained at the Keck telescope in July (H) and September (J and Ks) 2008. 
The inner part of the H-band image has been rotated by 1 degree to compensate for the orbital 
motion of the d between July and September. The central region is masked out in the upper 
images but left unmasked in the lower to clearly show the speckle noise level near d. 
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Figure 2: HR 8799bcd astrometric analysis. The positions of HR 8799bcd at each epoch are 
shown in both the overall field of view and in the zoomed-in insets. The solid oscillating line 
originating from the first detected epoch of each planet is the expected motion of a unbound 
background objects relative to the star over a duration equal to the maximum interval over which 
the companions were detected (4 years for b and c, two months for d.) All three companions 
are confirmed as co-moving with HR 8799 to 98cr for b, 90a for c and ~ 6a for d. Counter- 
clockwise orbital motion is observed for all three companions. The dashed lines in the small 
insets connect the position of the planet at each epoch with the star. A schematic dust disk - 
at 87 AU separation to be in 3:2 resonance with b while also entirely consistent with the far- 
infrared dust spectrum - is also shown. The inner gray ellipses are the outer Jovian-mass planets 
of our Solar system (Jupiter, Saturn, Uranus & Neptune) and Pluto shown to scale. 



18 



10 


1 1 1 1 1 1 1 

- □ HR8799b 
OHR8799C 
O HR8799d 

; (field) 

+ Pleiades 


1 1 1 1 1 1 1 1 ii 1 1 1 1 1 1 1 1 1 1 1 1 

▲ 

AB Pic b \ t 


12 




• 


14 
16 


• 


• 

• # ' 5 ' 2M1207 b - 
• • 

• ;• 

• • 




• 


• 


18 


1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



-1.0 -0.5 0.0 0.5 1.0 1.5 

H - K 

Figure 3: Absolute magnitude in H-band versus H-K color. Old, field (grey dots) and young 
Pleiades brown dwarfs (plusses) are shown along with 2 very low-mass brown dwarfs/planetary 
mass companions (filled black symbols). Open symbols are HR 8799b (square), c (diamond), 
and d (circle). 
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log (age) [yr] 

Figure 4: Luminosity versus time for a variety of masses [34]. The three coeval points are 
HR 8799b (square), c (diamond), and d (circle); c and d data points are displaced horizontally 
for clarity. The locations of the low mass object AB Pic b on the planet/brown dwarf dividing 
line and a planetary mass companion (2M 1207b) to the brown dwarf 2M1207 are also shown 
(note that alternative models proposed for 2M1207 lead to somewhat larger luminosity and mass 
(~8 Mjup) for the companion p2l). The deuterium burning mass limit, currently believed to be 
~13.6 Mjup, has been incorporated into a "working definition" of a planet by the International 
Astronomical Union and is used here to separate planets (which also must orbit a star) from 
brown dwarfs. The boundary between stars and brown dwarfs is set by stable hydrogen burning. 
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1 2 3 4 



Wavelength (/xm) 

Figure 5: Synthetic spectra from model atmospheres containing clouds located between 10 and 
0. 1 bar of pressure are compared to the measured fluxes (with 3 sigma error bars) for HR 8799 b, 
c and d. Response curves for each filter band pass are indicated along the x-axis. The predicted 
magnitudes from the synthetic spectra, averaged over the filter passbands, are shown by the 
filled symbols. 
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Supplemental Online Material 
ADI Observations and Data Reduction 

The first ADI imaging sequence on HR 8799 was obtained on October 17th, 2007 at the Gem- 
ini Observatory (program GN-2007B-Q-77) using the ALTAIR AO system (SI) and the NIRI 
camera (S2). After the first set of observations, follow-up multi-band imaging was carried out 
at the Gemini telescope (program GN-2008B-Q-64) and the W.M. Keck II telescope with its 
AO system (S3) along with the NIRC2 narrow-field camera (S4). To derive accurate bolometric 
luminosities, we have acquired a set of J-, H-, K'-, Ks- and L'-band images that span most of 
the spectral energy output of the companions. In addition, we have acquired a set of CH4S- 
and CH4L-band data (two filters that are in and out of the 1.6/xm methane absorption feature 
seen in cloud-free cool atmospheres) to search for a methane spectral signature. All the data 
were acquired in the linear regime (< 4%) of the detector. Once the planetary companions were 
noted in the 2007-2008 data we re-analyzed non-ADI NIRC2 data obtained on July 10 2004 at 
the Keck Observatory for a related project, and detected the two outermost companions of the 
system. In all data sets, both the 38 and 68 AU companions are detected. The 24 AU companion 
is only detected in the most sensitive ADI image sets. 

The data reduction for all ADI sequences are performed using a custom DDL script following 
a standard reduction technique (55). A shutter-closed dark image having the same observing 
parameters is first subtracted to remove the detector electronic bias and the image is divided by 
a flat field (produced by uniform illumination) to normalize the spatial variations in sensitivity. 
For the longer- wavelength K', Ks and L' images, we subtract a blank-sky image to remove 
thermal backgrounds. Deviant pixels are then removed by interpolating adjacent pixels and 
image distortions are corrected (for Gemini, the IDL observatory script is used, while for Keck 
the IDL script nirclwarp.pro is applied). We coarsely register the images using the unsaturated 
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data and a cross-correlation technique is used for fine registration of all the saturated/occulted 
images to sub-pixel accuracy. A 21 x 21 pixel unsharp mask filter is then applied to remove 
the low-spatial frequency noise (this filter is applied to both the unsaturated and saturated data). 
The improved ADI LOCI data reduction scheme (S6) is used for each individual image to select 
a set of reference images in which the stellar PSF is sufficiently similar but the field-of-view 
has rotated enough to not attenuate the companions. After subtracting the reference from each 
frame we de-rotate and average the frames to produce the final combined image. Fig. 1 shows 
the resulting data reduction for a subset of the HR 8799 acquired data. 

The exoplanet relative positions are obtained at each epoch by fitting a Gaussian function to 
the planet intensity profile. If multiple data sets are available for a given epoch, the position of 
the exoplanets are averaged together to reduce random errors. The measured positions were then 
compared with the expected drift for an unbound background/foreground object from HR8799's 
space motion and the Earth-induced parallax effect. Due to the large number of data sets avail- 
able, the higher signal-to-noise ratio of individual sets and well characterized astrometry, only 
the Keck data has been considered to perform the proper motion analysis. The relative detector 
positions of the exoplanets and the star are transformed to sky coordinates using the NIRC2 
narrow camera plate scale value (9.963 ± 0.005 mas/pixel) and North orientation [0.13 ± 0.02 
degree, S7]. Due to the uncertainty in registering the saturated images at the image center to 
perform the ADI processing and field rotation alignment, a conservative 0.5 pixel (~5 mas) 1 
sigma centroid uncertainty is added for each data set at each epoch. Note that uncertainty in 
both the plate scale and North orientation are well below the calculated centroid accuracy of 
our companions. Since all astrometric observations were acquired using NIRC2, whose plate 
scale is known to be very stable with time [<0.02% variations of plate scale since 2004, 57], 
the absolute plate scale and orientation errors are not included in the positional uncertainties 
given here. If the Keck astrometry is compared to that from other telescopes, those uncertain- 
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ties would need to be included. Table S2 summarizes HR 8799bcd astrometry for all acquired 
data sets. 

Accurate exoplanet photometry is essential to derive the companions' basic physical charac- 
teristics. Due again to higher SNRs and the fact that d is detected at almost all wavelenghs, the 
Keck data is preferred for the photometry analysis. For J-, Ks- and L'-band, we have chosen the 
coronagraphic data over the earlier saturated PSF sequences due to the better photometric cali- 
bration that it provides. The exoplanets' fluxes are derived after convolving the combined ADI- 
reduced images by a circular aperture having a diameter equal to the wavelength-dependent 
fuU-width-at-half -maximum of the telescope diffraction limit. The relative star-to-planet rela- 
tive intensity is obtained by taking the ratio of the peak fluxes (after aperture convolution) of 
the average unsaturated stellar images to the planetary fluxes. For coronagraphic long exposure 
data, the relative planet-star photometry is obtained from the unsaturated PSF core detected 
through the partially transmissive focal plane mask. The coronagraph throughput has been cal- 
ibrated using unsaturated unocculted images. The ADI processing does somewhat attenuate 
the fluxes from the companions. The stellar PSF obtained from the unsaturated data was used 
to introduce artificial sources into the raw images at various separations, position angles and at 
intensities similar to the ones of the detected exoplanets. These artificial images were processed 
using the same pipeline to calibrate the ADI algorithm throughput as a function of position (the 
Fig. 1 detection image has been renormalized to show a throughput of 1 at all separations). In- 
dividual magnitudes for each companion were derived by calculating the relative intensity (for 
each bandpass) of the planets compared to that of the host star. The star infrared magnitudes 
were taken from 2MASS and its L' magnitude was derived from S8. The derived photome- 
try (absolute magnitudes) at all epochs and wavelengths can be found in Table 1. Photometric 
errors are estimated from several sources. First, the variations in peak intensity of successive 
images of the unsaturated PSF (typically ~ 5%). Second, the speckle and photon noise in the 
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scattered starlight halo, estimated from image statistics at the separation of each companions 
(data set and separation dependent, from 1% to 43%). Third, the ADI-induced photometric er- 
ror from possible registration offsets from the exact rotational center of the field of view (~ 0.5 
pixel RMS resulting in a 3% error). Fourth, the residual registration accuracy between images 
(~ 0.5 pixels inducing a 3% error) and finally the ADI flux renormalization error (typically 
~ 5%) estimated from the artificial sources. For non-coronagraphic ADI sequences, an addi- 
tional 15% photometric error is added to account for potential seeing variations during long 
saturated image sequences (this fluctuation was estimated from coronagraphic data). 

Planet Atmosphere Modeling 

To determine the luminosity of each planet, synthetic spectra from a variety of planetary atmo- 
sphere models were fit to the photometric measurements described above. These models and 
spectra were generated using the PHOENIX atmosphere code (S9) which includes two methods 
for describing the impact of cloud formation from solid (and liquid) material that condenses 
under the assumption of chemical equilibrium. These two cloud models cover the cases of 
extremely cloudy, with nearly uniform coverage both vertically and horizontally in the atmo- 
spheres, and completely cloud-free (SIO). It is known from studies of brown dwarfs (which 
experience similar atmospheric conditions as giant planets) that atmospheres transition from 
having very cloudy atmospheres to nearly cloud-free atmospheres as they cool (called the L-T 
transition). Consequently, cloud models that are intermediate between the two extreme cases 
are required to reproduce the diversity observed in ultra-cool objects. Since the planets orbiting 
HR 8799 have photometric properties similar to brown dwarfs at or near the L-T transition, nei- 
ther of the two extreme cloud models are likely to be appropriate. To better describe the planets, 
an intermediate cloud model similar to that described by Sll was added to PHOENIX. In this 
case, the clouds are more confined vertically in the atmosphere than in the extreme DUSTY 
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case, where the number density of cloud particles decreases with increasing height. While this 
cloud model is certainly more realistic than the extreme cloudy assumption, modeling clouds in 
giant planet atmospheres is very difficult and a wide variety of theoretical attempts are actively 
being explored and, so far, no clear physical picture has emerged (S12). 

The synthetic spectra from these models were convolved with the response curves of each 
near-IR filter (Fig. 5, lower panel) used for the observations, thereby generating synthetic pho- 
tometry. A standard Levenberg-Marquardt least-squares minimization procedure was used to 
fit these synthetic photometry to the real data. As already mentioned in the main article, the 
best fit temperatures and radii for cloud assumptions were very different, while the luminosities 
were very similar. This is to be expected since the observations cover a large fraction of the 
total spectral energy distribution. As a further check on the luminosity, black bodies and semi- 
empirical bolometric corrections were used and resulted in nearly identical values (see Table 1 
for the luminosities and uncertainties). 

Interestingly, the best fit effective temperatures from the model atmosphere analysis are 
high (ranging from 1700 to 1400K). Such high effective temperatures would imply either un- 
realistically small radii (a few tenths the size of Jupiter) or that each planet has its own disk 
exactly aligned edge-on with our line-of-sight causing several magnitudes of extinction. This 
later situation is very unlikely since edge-on disks are rare, and the observed orbital motions 
along with the small projected rotation of the host star suggest that HR 8799 is viewed closer to 
face-on. The most likely explanation for the unreasonably high effective temperatures implied 
by the best-fit atmosphere models is missing atmospheric physics (e.g., poor cloud modeling 
and possibly non-equilibrium chemistry). 

Since the theoretical cooling tracks are far less sensitive to the assumptions regarding at- 
mospheric clouds, the masses, effective temperatures, and radii for each planet were obtained 
from cooling tracks (see Fig. 4) rather than the best fit values to the synthetic photometry men- 
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tioned above. However, as a consistency check, synthetic planet spectra were computed using 
our intermediate cloud model and the parameters determined from the cooling tracks. These 
synthetic spectra are compared to the observed data in Fig. 5 and while they do not represent a 
best-fit to the data, they are a reasonably close match and lend confidence to the values implied 
by the luminosity and cooling track comparison. Some of the discrepancies with these synthetic 
spectra could be alleviated by including non-equilibrium chemistry which has been shown to 
impact the strength of CH4 and CO bands (S13). 
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SOM Tables 



Table SI: HR 8799 observing log 



Inst./Telescope 


UT Date 


Filter 


Saturated Data 
Total Exp. Time (s) 


NIRC2/Keck 


2004 July 14 


H 


480 




2007 Oct. 25 


CH4S 


2700 




2008 July 11 


H 


1740 






CH4S 


1680 




2008 Aug. 12 


J 


900 






H 


420 






CH4L 


870 






K' 


840 






L' 


900 




2008 Sept. 18 


J 


1370 






Ks 


1200 






U 


2800 


NIRI/Gemini 


2007 Oct. 17 


CH4S 


3600 
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Table S2: NIRC2 HR 8799bcd astrometry. 



UT Date Filter Separation w.r.t the host star in [E, N] 







b 










2UU4 July 14 


TT 
M 


[i.4 < i, U.oo4J 


F n VQO n fii ol 
[— U. / cjy, U.oizJ 








zUU/ UCt. ZD 




[i.oiz, U.oUoJ 


[— U.D / 4, U.DoiJ 








zooo July 1 1 


TT 

M 


[i.ozo, U.oUoJ 


[— U.DiDD, U. / UoJ 


r 

[- 


-U.ZU4, 


— U.o /yj 






[l.i3Z0, U. /yzj 


F n fifin n find 
[— u.oou, u.oyoj 


F 

[- 


-U.zii, 


— U.oooJ 


Avg. (±U.UU4 ) 




[i.oz 1, u. / yyj 


F n fi n vni 1 
[— U.OOo, U. / UIJ 


r 

[- 


-U.zUo, 


— U.oozj 


ZUUO Aug. iZ 


T 
J 


Fi c;9fi n sni 1 

[i.OZO, U.oUiJ 


[ — u.ooy, u. ( uzj 










H 


[1.527,0.803] 


[-0.657,0.708] 










CH4L 


[1.522,0.797] 


[-0.662, 0.703] 


[- 


-0.217, 


-0.582] 




K' 


[1.527,0.802] 


[-0.657,0.703] 


[- 


-0.213, 


-0.580] 




L' 


[1.531,0.802] 


[-0.652,0.713] 


[- 


-0.217, 


-0.584] 


Avg. (±0.002") 




[1.527,0.801] 


[-0.657,0.706] 


[- 


-0.216, 


-0.582] 


2008 Sept. 18 


J 


[1.527,0.796] 


[-0.656,0.705] 


[- 


-0.213, 


-0.582] 




Ks 


[1.525,0.796] 


[-0.659, 0.708] 


[- 


-0.221, 


-0.581] 




L' 


[1.532,0.802] 


[-0.657,0.703] 


[- 


-0.214, 


-0.584] 


Avg. (±0.003") 




[1.528,0.798] 


[-0.657, 0.706] 


[- 


-0.216, 


-0.582] 
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